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ABSTRACT 

The  effect  on  cable  reliability  of  random  cyclic  loading  such  as 
that  generated  by  the  wave  induced  rocking  of  ocean  vessels  deploying 
these  cables  is  examined.  A simple  model  yielding  explicit  formulas  is 
first  explored.  In  this  model,  the  failure  time  of  a single  element  under 
a constant  load  is  assumed  to  be  exponentially  distributed,  and  the  random 
loadings  are  a two  state  stationary  Markov  process.  The  effect  of  load 
on  failure  time  is  assumed  to  follow  a power  law  breakdown  rule.  In  this 
setting,  explicit  results  concerning  the  distribution  of  bundle  or  cable 
failure  time,  and  especially  the  mean  failure  time,  are  obtained.  Where 
the  fluctuations  in  load  are  frequent  relative  to  cable  life , such  as  may 
occur  in  long-lived  cables , it  is  shown  that  randomness  in  load  tends  to 
decrease  mean  cable  life , but  it  is  suggested  that  the  reduction  in  mean 
life  often  can  be  restored  by  modestly  reducing  the  base  load  on  the  structure 
or  by  modestly  increasing  the  number  of  elements  in  the  cable. 

In  later  pages  this  simple  model  is  extended  to  cover  a broader  range 
of  materials  and  random  loadings.  Asymptotic  distributions  and  mean  failure 
times  are  given  for  cable  elements  that  follow  a Weibull  distribution  of 
failure  time  under  constant  load,  and  loads  that  are  general  nonnegative 
stationary  processes  subject  only  to  some  mild  condition  of  asymptotic 
independence . When  the  power  law  breakdown  exponent  is  large , the  mean 
time  to  cable  failure  depends  heavily  on  the  exact  form  of  the  marginal 
probability  distribution  for  the  random  load  process  and  cannot  be  summarized 
by  the  first  two  moments  of  this  distribution  alone. 
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The  Time  to  Failure  of  Cables  Subjected  to  Random  Loads 


1.  Introduction. 

Wl/VWVVV\AVWV\/ 

The  Naval  Undersea  Center  in  Hawaii  is  involved  in  the  construction  of 
electromechanical  cables  for  deep  sea  operation,  some  of  which  are  20,000 
feet  long.  Kevlar-49^  is  being  used  in  the  strength  members  because  of  the 
tremendous  weight  savings  over  steel,  which  can  hardly  support  its  own 
weight  in  these  lengths.  The  severing  of  such  cables  is  enormously  expensive 
because  of  the  equipment  involved  on  the  ocean  floor. 

This  paper  examines  the  effect  on  cable  reliability  of  random  cyclic 
loading  such  as  that  generated  by  the  wave  induced  rocking  of  ocean  vessels 
deploying  these  cables.  A simple  model  yielding  explicit  formulas  is  first 
explored.  In  this  model,  the  failure  time  of  a single  element  under  a 
constant  load  is  assumed  to  be  exponentially  distributed,  and  the  random 
loadings  are  a two  state  stationary  Markov  process.  The  effect  of  load  on 
failure  time  is  assumed  to  follow  a power  law  breakdown  rule.  In  this 
setting,  explicit  results  concerning  the  distribution  of  bundle  or  cable 
failure  time,  and  especially  the  mean  failure  time,  are  obtained.  Where 
the  fluctuations  in  load  are  frequent  relative  to  cable  life,  such  as  may 
occur  in  long-lived  cables,  it  is  shown  that  randomness  in  load  tends  to 
decrease  mean  cable  life,  but  it  is  suggested  that  the  reduction  in  mean 
life  often  can  be  restored  by  modestly  reducing  the  base  load  on  the  struc- 
ture or  by  modestly  increasing  the  number  of  elements  in  the  cable. 

In  later  pages  this  simple  model  is  extended  to  cover  a broader  range 
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of  materials  and  random  loadings.  Asymptotic  distributions  and  mean  failure 
times  are  given  for  cable  elements  that  follow  a Weibull  distribution  of 
failure  time  under  constant  load,  and  loads  that  are  general  nonnegative 
stationary  processes  subject  only  to  some  mild  condition  of  asymptotic 
independence.  Cable  behavior  in  these  more  general  circumstances  may 
differ  from  that  found  in  the  simple  model  first  explored.  The  exact  form 
of  the  marginal  distribution  of  the  stationary  load  process  seems  critical 
when  the  exponent  in  the  power  law  breakdown  rule  is  large. 
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2.  A simple  model  with  deterministic  loads. 

V\AAAAA/w'yWV\AA>VWWWWWW\/VWWWWVW\/ 

A single  fiber  subjected  to  the  time  varying  tensile  load  £(t)  fails 
at  the  random  time  T.  We  postulate  the  failure  time  distribution 

Pr{T  < t)  = 1 - exp{-  K[i,(s)]ds), 

which  corresponds  to  the  failure  rate  or  hazard  rate  r(t)  = K[Jl(t)]. 

That  is,  a single  fiber,  having  not  failed  prior  to  time  t and  carrying 
load  £(t),  will  fail  during  the  interval  [t,  t+At)  with  probability 
K[i.(t)]At  + o(At)  where  o(At)  denotes  remainder  terms  of  order  less  than 
At  as  At  vanishes. 

The  function  K,  called  the  breakdown  rule,  expresses  how  changes  in 
the  load  affect  the  failure  probability.  We  concentrate  on  the  power  law 
breakdown  rule  in  which  K( £ ) = k£P  for  positive  constants  x and  p. 
Under  a constant  load  d(t)  = l,  the  failure  time  of  a fiber  obeying  power 
law  breakdown  is  exponentially  distributed  with  mean  E[ T | £.3  = 1/K(£)  = 
£-p/x.  A plot  of  mean  failure  time  versus  load  is  linear  on  log-log  axes, 
a relationship  which  is  commonly  observed  in  fatigue  and  stress  rupture 
studies. 

Phoenix  (1976)  gives  the  following  estimates  for  p: 


Strand  type  p_ 
Kevlar-49/Epoxy  42 
Graphite  Fiber/Epoxy  78 
S- Glass/Epoxy  30 
Beryllium  Wire/Epoxy 
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The  analysis  of  cable  behavior  when  p < 1 involves  several  additional 
technical  nuisances  that  do  not  significantly  add  to  our  understanding  and 
that  are  not  present  when  p >.  1.  In-as-much  as  Phoenix's  estimates  indi- 
cate that  a restriction  to  p ^ 1 is  not  severe  in  practice  for  some  material 
types,  henceforth  we  consider  only  values  p >.  1. 

Now  place  n of  these  fibers  in  parallel  and  subject  the  resulting 
bundle  or  cable  to  a total  load,  constant  in  time,  of  nL  units,  where  L 
is  the  nominal  load  per  fiber.  What  is  the  probability  distribution  of  the 
time  at  which  the  bundle  fails?  In  the  next  section  we  will  see  that  the 
affect  of  the  constant  load  L and  the  power  law  breakdown  coefficient  k 
can  both  be  absorbed  in  a simple  scale  change.  Therefore  in  the  remainder 
of  this  section  we  concentrate  on  a unit  load  L = 1 and  k = 1. 

Since  the  fibers  are  in  parallel,  the  bundle  failure  time  equals  the 
failure  time  of  the  last  fiber.  Let  be  the  time  that  the  first 

(earliest)  fiber  fails,  let  Sj  be  the  time  that  the  second  fiber  fails, 
and  so  on.  The  bundle  carries  total  load  nL  = n,  and  at  the  start,  assuming 
equal  load  sharing,  each  of  the  n fibers  carries  load  nL/n  = 1.  Any 
particular  fiber  has  the  survival  distribution 

Pr(T  > t>  = e“K(1)t;  t >0, 

and  S^,  being  the  minimum  of  n such  individual  fiber  failure  times,  has 
a survival  distribution  which  is  the  nth  power  of  this, 

Pr{Sx  > t>  = e"nK(1)t,  t > 0. 

That  is,  has  an  exponential  distribution  with  parameter  nK(l). 
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When  the  first  fiber  fails,  each  of  the  n-1  remaining  fibers  carries 
load  nL/(n-l)  = n/(n-l)  and  has  the  failure  rate  K[n/(n-l)J.  Since  there 
are  n-1  fibers  remaining  in  the  bundle,  the  rate  at  which  the  next 
failure  in  the  bundle  occurs  is  (n-1)  times  this,  or  (n-l)K[n/(n-l)], 
which  leads  to 


MS2  - sx  » t> 


-(n-l)K[n/(n-l)Jt 


This  reasoning  continues,  and  we  deduce  that,  when  i fibers  have  failed, 
the  n-i  remaining  fibers  each  carry  load  n/(n-i)  and  have  individual 
failure  rates  K[n/(n-i)].  The  rate  at  which  the  next  failure  occurs  among 
the  n-i  fibers  remaining  in  the  bundle  is 


Pr(S.  , - S.  > t}  = e-(n-i)KCn/(n-i)]t 
l+l  1 


t > 0. 


This  analysis  allows  us  to  write  the  bundle  failure  time  as  a sum 

of  independent  and  exponentially  distributed  differences  = sn  " l’ 

Y2  = Sn-1  ' Sn-2”--’Yn  = S1  " S0  where  S0  = °*  That  is 


where 


S_  = Y, 


+ Y 

n 


Pr{Yk  > t}  = exp{-  nPk1_Pt),  t >,  0. 


We  have 


L 
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ECYk]  = n~pkP_1, 

Var[Yk]  = n"2p  k2p_2, 

ECExp{-sYk}]  = npk1'p/(s+npk1"p) 

= 1/(1  + skP  E/nP ) for  s > -npk1  P, 


and 

ECSn]  = ^k=l  <k/n>P"1<l/n>  (2.D 

Var[S  ] = (1/n)  (k/n)2p_2 (1/n ) (2.2) 

n K-i 

E[exp(-sS  }]  = 1/nP.Hl  + s(k/n)P'1(l/n)],for  s > -1.  (2.3) 

n K- j. 

Asymptotic  distribution.  Study  of  the  limiting  distribution  of  as  n 

becomes  large  divides  into  several  cases  according  to  the  value  of  p.  The 
simplest  case  is,  at  the  same  time,  the  most  important  case  in  practice  and 
concerns  p >.  1.  We  concentrate  on  it.  Then 


H«  KSn]  - Urn  I"=1  ± 

n-*»  n-x» 

= xP_1  dx  = 1/p, 

lim  n Var[Sn]  = lim  j£=1  ^)2f>  2 

n-x»  n-x» 

= /J  X2p_2  dx  = 1/ ( 2p-l ) , 

which  suggests  that  2^  = *^n(Sn  - 1/p)  should  have  a limiting  distribution 
whose  mean  is  zero  and  whose  variance  equals  1/(2 P-1).  Indeed,  this  is  the 
case,  the  limiting  distribution  being  normal  (Gaussian).  This  may  be 


established  via  Liapounov's  central  limit  theorem  or  by  a direct  argument 
invoking  Levy’s  convergence  theorem  for  moment  generating  functions.  We 
proceed  with  the  latter  approach.  First 
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E[exp(-s/n(S  - £ (^)P  1 = 

k=l 

exp(s/n  (— )P  1 — } 
k=l  r n n 

n"  (1  + s/n  (-)p_1 
k=l  n n 


TT{- 

k=l  v>- 


k.p- 


+ s/n  (— ) 
n 


11  12  ,k . 2p-2  1 -3/2. 

— + tts  (— ) — + 0(n  ) 

n 2 n n 


1 + s/n  (— )P  1 — 


= YJ  u + fs2(|)2p'2  i + 0(n"3/2)} 

k=l 

— » expt-^s2  x2p  2 dx)  = exp{-|s2/(2p-l) }. 

n-H»  r 2 0 i 

The  limit  above  is  an  easy  consequence  of  the  lemma  preceding  Theorem  7.1.2 

in  Chung  (1968).  As  the  limit  function  is  the  moment  generating  function 

corresponding  to  a mean  zero  normal  distribution  having  variance  l/(2p-l), 

this  must  be  the  asymptotic  distribution  for  /n  (S  -7"  (— )p_1  -i) . But 

n ‘■In  n 

\ - ^ \ 

where 




R = /nil1}  (-)p_1  - - 1/p), 
n ^1  n n 

r“/rn  /^nP-1  1 rl  P-1  j n 
= ^(Il  (n}  n"/oX  dx)- 
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, , .k.p-l  p-1  ,k-l,p-l  k k-1 

Because  p > 1 we  have  (— )K  > x > ( y whenever  — > x > , 

= n = — n n—  — n 

whence 

0 < R < & l {(-)P_1  - - (— — -)P_d  -} 
k=l 

= /n  (^-)  (Collapsing  sum.) 


Thus  R -*■  0 as  n -*■  00  and  it  follows  by  Slutsky's  theorem  (Cramer,  1945 
n 

Theorem  20.6,  p.  254)  that  ZR  = /n  (S^  - 1/p)  asymptotically  shares  the 
same  distribution  as  in  (Sn  - (^-)P  1 ^-) . That  is,  is  asymptotically 

normally  distributed  with  mean  zero  and  variance  l/(2p-l)  as  claimed. 


It  is  often  suggestive  to  write  this  result  in  the  form 


S = 1/p  + Z /in 
n n 


where  Zfi  has  a limiting  zero  mean  normal  distribution  with  variance  equal 
to  l/(2p-l) . 

For  what  follows,  let  us  define  M(s)  to  be  the  number  of  unfailed 
fibers  at  time  s assuming  the  unit  nominal  load  per  fiber  L = 1 and 
power  law  coefficient  k = 1 . Then 


M(s)  = n for  0 4 s < S^, 


= n-1  for  < s < Sj, 


=1  for  S < s < S , 

n-1  — n’ 

= 0 for  S 4 s, 

n 


and  Y is  the  sojourn  time  for  the  process  in  state  k.  The  independent 
k 

exponential  distributions  for  Y^,  Yn_^,...,Y^  show  that  (M(s);  s ^0} 
is  a pure  death  stochastic  process  with  death  rates 

Pr{M(t+At ) = k-1 ] M(t ) = k}  = npki  pAt  + o(At), 

and  S = inf{s  >0;M(s)=0}  is  the  first  time  this  process  hits  zero, 
n — 

Time  varying  loads.  Consider  now  a bundle  of  n fibers  following  power 
law  breakdown  K(K,)  = kIp  and  subjected  to  a total  load  nL(t)  that  may 
vary  with  time.  To  avoid  trivialities,  assume  that  L(t)  is  strictly 
positive  for  positive  t,  define 

H(w | L)  = kJq  L(t)PdT 


and  let 


G(sjL)  = H 1(s | L) 

be  the  inverse  to  the  strictly  increasing  function  H( » | L) . For  typographi- 
cal convenience  we  will  often  omit  the  load  L from  our  notation  and 
simply  write  H(w)  for  H(w)l)  and  G(s)  for  G(s|L). 

We  will  relate  the  time  varying  load  problem  to  the  constant  load 
problem  by  using  H as  a time  scale  change.  Let  N(t)  be  the  number  of 
unfailed  fibers  at  time  t.  Introduce  the  rescaled  process 


M(s)  = N(t)  where  s = H(t|L). 
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Evaluate 

Pr{M(s+ds)  = k-1 | M( s ) = k}  = Pr{N(t+dt)  = k-l|N(t)  = k} 

= k k (nL(t)/k)P  dt 
= r^k1  pds,  since  ds  = <L(t)pdt, 

to  discover  that  M(s)  evolves  as  a bundle  of  fibers  carrying  unit  nominal 
load  and  for  which  k = 1.  If  Wn  denotes  the  bundle  failure  time  under 
the  time  varying  load,  then  Sn  = H(Wn|L)  is  the  failure  time  of  a bundle 
subjected  to  a constant  unit  load  and  we  have  the  explicit  respresentaticn 

W = g(sJl) 

n n ' 

which  relates  the  bundle  failure  time  Wn  under  varying  load  and  arbitrary 
k to  the  failure  time  under  unit  load  and  with  < = 1. 

Under  a constant  nominal  load  L(t)  H L we  have  H(w)  = <LPw,  and 
G(s)  = s/(<L  ) so  that  = S^/ ( kLP ) and 

W = l/(p<LP)  + Z,/*/n,  (2.4) 

n n 

where  Z*  = Z /(kLP)  is  asymptotically  normally  distributed  with  mean  zero 
n n 

and  variance  equal  to  1/[kLP  /(2p-l)]^. 

This  analysis  of  bundle  failure  time  under  deterministic  loads  and 
the  assumed  exponentially  distributed  failure  times  of  single  fibers  under 
constant  load  is  due  to  Coleman  (1958),  and  later,  to  Birnbaum  and  Saunders 
(1958). 


11 


While  explicit  results  concerning  bundles  of  finite  size  n are  not 
generally  obtainable,  the  asymptotic  normal  distribution  of  bundle  failure 
times  prevails  under  more  general  circumstances  than  so  far  indicated.  In 
particular,  suppose  that  a single  fiber  subjected  to  a load  £(s)  for 
s ^ 0 follows  the  failure  time  distribution 

F(t  1 i,)  = 'K/J  <8.(s  )pds ) 

for  some  cumulative  distribution  function  y.  Introduce  the  notation 

g(u)  = y "^u)  for  the  inverse  function  to  y,  and  <f(y)  = (l-y)p.  Next, 

let  S be  the  failure  time  of  a bundle  of  such  fibers  subjected  to  a unit 
n 

nominal  load  per  fiber  and  for  which  k = 1.  For  a broad  class  of  functions 
y,  Phoenix  (1977)  has  shown  that,  asymptotically  for  large  n,  the  standard- 
ized bundle  failure  time  2 ; = >/n{Sn-y}  is  normally  distributed  with  mean 

zero  and  variance  o2 , where 

p = Jq  4>(y)g'  (y)dy 

o2  = Jq  (u)<f>'  (v)g'  (u)g'  (v)(uAv-uv)dudv, 

with  "prime"  denoting  differentiation. 

The  most  interesting  special  case  concerns  the  particular  function 
y(x)  = 1 - exp{-cxu)  where  c and  a are  fixed  positive  parameters. 

In  this  case,  the  failure  time  T of  a single  fiber  subjected  to  a constant 
load  £(t)  = L follows  the  Weibull  distribution 


F(t  |l)  = 1 - exp(-  c(KLPt)a),  t > 0 
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and  then 


a2  = 2( 


1/a  p 
ac  kL 


1 rl,,  -^p-1  r , 1 , , (l-a)/a , 

^ = “17^  Jo(1"y)  u°s(ir7)}  dy’ 

ac 

>2  /Jd-v)|,-1uog(Tir))a-“)/7jJa-u)('-2(iog(IiI)}<1-“)/«<1u  dv. 


The  Weibull  family  possesses  several  properties  that  are  desirable  for 

distributions  of  time  to  failure  or  strength  of  filaments.  In  particular, 

if  the  time  to  failure  of  a typical  segment  of  a filament  follows  a Weibull 

distribution,  then  the  time  to  failure  of  a filament  comprised  of  several 

independent  segments  will  also  follow  a (different)  Weibull  distribution. 

For  a review  of  the  Weibull  distribution  in  this  context,  see  Harlow, 

Smith  and  Taylor  (1978),  which  contains  other  material  as  well. 

When  c = a = 1 we  return  to  the  exponential  case  treated  in  such 

detail  in  the  beginning  of  this  section,  and  then  p = 1/p  and 
2 

o = l/(2p-l),  as  may  be  verified  easily. 

Returning  to  the  general  Weibull  distribution  ¥(x)  = 1 - exp{-cxa), 
if  W^  is  the  failure  time  of  a bundle  having  the  time  varying  nominal 
load  per  fiber  of  L(s),  s >_  0,  and  whose  breakdown  rule  has  an  arbitrary 
value  k , then  again  we  have  the  representation 


W = G(S  |L) 
n n 1 


The  crucial  assumption  here  is  the  postulated  power  law  breakdown  rule 


K(fc)  = <£P. 
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Given  L(0)  = L+A,  the  process  H(w)  increases  linearly  at  rate  a+ 
for  a duration  which  is  exponentially  distributed  with  parameter  X.  It 
then  increases  linearly  at  rate  a_  for  a similarly  distributed  time,  and 
so  on.  The  inverse  function  G(s)  = H 1(s)  has  a similar,  but  distinct, 
behavior.  It  increases  at  rate  l/a+  for  an  exponentially  distributed 
time  having  parameter  A/a+.  It  then  increases  at  rate  l/a_  for  an 
exponentially  distributed  time  having  parameter  X/a_  , and  so  on.  We  have 
the  representation 


G(s)  = fSQ  Y(o)do 


where  {Y(o);  o >_  0}  in  a two-state  Markov  process,  the  states  being 
l/a+  = 1/[k(L+A)p]  and  l/a_  = 1/[<(L-A)P]  and  the  times  in  these  states 
being  exponentially  distributed  with  respective  parameters  A/a+  and 
A/a_ . The  transition  function  for  Y is  known  to  be 


a_  a 1 i 

P (a)  = + — - — exp{-A(—  + — )a), 

a++a_  a++a_  r a+  a_ 


and 


a - a_  . , 

exp{-A(—  + — )o). 
_ _ a a 

+ - + - + 


P (a)  = — — 

+-  a +a  a +a 


(See  Karlin  and  Taylor  (1975)  Problem  7,  page  154.)  Since  Y(0)  = 
or  l/a_  with  probability  1/2  each,  we  have 


".PlWy* 


1 -w-> 


P (o)  = Pr(Y(a)  = 1/a  } 


= -±(P+_(a)  + P__(a) ) 
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a+  + a 


1 , + 
+ o (— 


2 a+  + a. 


-)exp{-X( — + — )o} 
+ 


and 


P+(o)  = 1 - P_(o) 


a+  + a 


- i (*♦ 

n v 


2 'a+  + a“)exp{-A(i;  + ^)a)> 


Then 


E[Y(o)]  = — P (o)  + i-  p (o) 


a+  + 


1,1  1 a+  ‘ a 


T~r  + 2 (~  - r)(rrr)expt4(r  + s-)o}> 


+ + 


a+  a- 


. -CO 

which  we  abbreviate  in  the  form  E[Y(o)]  = A + Be  with 


A = 2/(a+  + a_), 

B = i (f  . f)(!i 


a+  + a 


(3.1) 

(3.2) 


and 


c = * r>- 


(3.3) 


Then,  since 


S 

Wn  = /0nY(o)do 


(3.4) 


we  have 


E[W  ] = E{E[WlS  ]} 
n n1  n 

S 

= E{/Qn(A  + Be-C°  )do} 

-CS 

= A E[S  ] + (B/CKl  - E[e  n]}.  (3.5) 

n 

In  conjunction  with  equations  (2.1)  and  (2.3),  this  provides  an  explicit 
formula  for  the  mean  time  to  bundle  failure  under  the  random  load  L. 

While  Equation  (3.5)  is  explicit,  it  never-the-less  remains  clumsy. 
The  formula  simplifies  in  certain  extreme  rases,  however,  which  enables 
us  to  gain  some  insight  as  to  how  random  loads  affect  mean  bundle  failure 
times.  When  X becomes  small,  the  frequency  of  the  load  fluctuations 
decreases  and  in  the  limit,  when  X = 0,  the  load  process  remains  at 
whichever  level  it  began.  We  then  have 


lim  E[W  ] 
XiO  n 


(A  + B )E[S  ] 
n 


= {- 


2 1,1 
+ a + 2 (a 


1 a+  ' a- 

- j-H—  ■--■~)}E[S  J. 
ci  ci  + 3 n 

+ + 


If  we  postulate  that  A is  small  compared  to  the  base  load  L,  and  that 
p is  moderate  so  that  0 < pA  <<  L,  we  have  the  Taylor  series  expansions 


a+  = tc(  L + A)P  = kLP(1  + pA/L  + P-(-P~^(^)2) 


and 
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lim  E[W  ] = ) { 

A + 0 kL°  1 + 


+ (T>2}E[sn] 


(-) 
2 vl' 


= (1/kLP  ) (1  + ip(p+l)(4/L)2}E[S  ],  0 < pA  <<  L.  (3.6) 


Now  it  is  doubtful  that  the  mean  failure  time  is  a relevant  criterion  of 
reliability  in  the  extreme  case  at  hand  since  half  the  time  the  bundle  will 
encounter  the  heavy  load  L + A and  will  fail  relatively  quickly.  Never- 
theless, recalling  that  the  mean  failure  time  under  no  random  load  (A  = 0) 
is  (1/kLP )E[Sn] , we  see  from  (3.6)  that  the  randomness  has  actually 
increased  the  mean  failure  time.  That  under  certain  circumstances  randomness 
could  actually  increase  mean  failure  times  was  a surprising  discovery. 

Much  more  relevant  in  practice,  it  would  seem,  is  the  case  where  X 
is  large  since  this  corresponds  to  many  random  fluctuations  in  the  load 
prior  to  bundle  failure.  Letting  A -►  °°  we  have 


lim  E[W  ] 
A-*»  n 


A E[Sn] 


7 E[SJ. 

a + a n 


When  p _>  1 we  have  the  convexity  inequality  that  (a+  + a_)/2  = 
(k/2 ) (L  + A)P  + {<17)  (L  - A)P  > <LP  whence 


lim  E[W  ] = - — l E[S  ] 

A—  n a+  + a.  n 

< (1/<LP)  E[S  ]. 
— n 


The  latter  being  the  strength  under  no  random  stress,  we  conclude  that 
rapidly  varying  random  loads  always  decrease  mean  failure  time. 
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If  we  introduce  the  notion  of  an  equivalent  load  L where 

eq 

rp  = — (L  + A)P  + 4 (L  - A)P  we  observe  that  a bundle  under  the  random 
eq  2 2 

load  L + A has  a mean  failure  time  equal  to  that  of  a bundle  under  the 

fixed  load  L „ when  A -*•  «. 
eq 

Let  us  highlight  from  (3.5)  that  the  design  formula  lim^^  E[W^] 

= A E[S  ] = l/(<pLP  ) is  conservative  in  the  sense  that  E[W  ] >_  A E[S  ] 
n eq  n — n 

for  all  A.  As  already  pointed  out,  the  presence  of  rapidly  fluctuating 
random  loads  tends  to  decrease  mean  failure  time,  and  it  is  of  major  interest 
to  estimate  the  additional  number  of  fibers  that  would  be  needed  to  restore 
the  loss.  An  increase  to  n^  fibers  from  nQ  fibers  will  decrease  the 
equivalent  load  to  (n(/ni^Leq  and  th:*-s  e<luated  to  the  deterministic 
design  load  L to  obtain 


nl  Leq 
n0  = L ' 


This  ratio  L /L  as  a function  of  A/L  for  various  values  of  p is 
eq 

evaluated  in  the  next  section,  along  with  the  same  ratio  under  some  different 
distributions  for  the  stationary  load  process. 
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It  is  possible,  on  the  one  hand,  to  extend  the  preceding  analysis  to 

cover  quite  general  stationary  load  processes,  and  on  the  other,  to  refine 

the  analysis  and  obtain  some  information  concerning  the  random  fluctuations 

in  failure  time  about  the  asymptotic  mean.  The  key  to  this  further  analysis 

is  a far  reaching  generalization  of  the  central  limit  theorem  due  to 

Stratonovich  (1968)  and  Has’minskii  (1966)  which  we  enunciate  only  in  a 

very  limited  case  that  suffices  to  meet  our  needs.  Let  (L(t);  t ,>  0}  be 

a stationary  process  taking  values  in  the  strictly  positive  bounded  interval 

Cl  . , L ].  Let 
mm  max 


L = E[L(t )]  and  L = {E[L(t)P]}1/p. 
eq 


Let  F^  be  the  o-algebra  of  events  generated  by  L(t)  for  s < t < t 
and  assume  an  asymptotic  independence  for  L(t)  in  the  form 


sup{  | P ( B ) A ) - P(B)|;  A 6 Fj,  B e F®^}  < B(s) 


where  8(s)  -*-0  as  s -*■  00  sufficiently  fast  so  that  s 8(s)  -»■  0 as  well. 

Let  F(£.)  be  a smooth  function  of  l,  bounded  on  the  interval  [L  . , L ] 

mm  max 

and  for  which  E[F(L(t))]  = 0,  and 


lim  | /T  J®  E[F(L(s))F(L(a))Jdo  ds  = 62  > 0. 
T-*® 


Next,  for  a fixed  A > 0,  consider  the  solution  Z^(t)  to  the  differential 
equation 
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dZ  (t)/dt  = /T  F(L(Xt ) ) , Zx(0)  = 0. 


Of  course,  the  solution  is  Z^(t)  = ^ Iq  F(L(Xs))ds  = (l/ZT)  / * F(L(v))dv. 

Note  that  large  values  for  X correspond  to  load  processes  L (s)  = L(Xs) 

having  "many  fluctuations"  prior  to  time  s,  or  alternatively,  observing 

the  fixed  load  process  L(v)  over  the  long  duration  Xt.  That  is,  large 

values  for  X are  quite  relevant  to  many  cable  design  situations.  The 

result  of  Stratonovich  and  Has'minskii,  restricted  to  the  case  at  hand, 

asserts  that,  as  X -*•  °°,  the  processes  (Z^Ct),  t >.  0}  converge  weakly  in 

the  space  of  continuous  functions  to  a Brownian  motion  process  having 

2 

variance  parameter  6 . 

We  apply  this  result  using  for  F the  smooth  function 


F(l)  = <(JlP  - LP  ),  L . <.  I < L 

eq  min  = — max 


whence 


Z, (t ) = L.(s)pds  - LP  t} 

A 4 0 A eq 


= /X  {H. (t ) - <LP  t} 
X eq 


where  Lx(s)  = L(Xs)  and  H^Ct)  = k /*  Lx(s)P  ds.  We  conclude  that 
{Z  (t),  t _>  0}  is,  asymptotically  for  large  X,  a Brownian  motion.  The 

A 

next  step  is  to  relate  this  convergence  to  the  behavior  of  Gx(s)  = H^Cs). 
To  this  end,  define 


V (s ) = /X  {G  (s ) - s/(kLp  )}. 
X X eq 
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i 

I 

J 


is  near  s/(k  Lp  ).  To  be  precise, 
eq 

bility  one  uniformly  on  compact  s 


these  convergences  take  place  with  proba- 
intervals.  We  exploit  this  by  writing 


V (s ) = -(1/k  LP  ) Z (s/<  LP  ) -(1/k  LP  ) R.(s) 
X eq  X eq  eq  X 


where 


R,  (s)  = Z.(G.(s))  - Z.  (s/k  Lp  ) 


X'  X 


eq 


We  claim  that  R (s)  converges  to  zero  in  probability  uniformly  on 

A 

compact  s intervals,  and  consequently,  that  V^(s)  converges  weakly  to 
a Brownian  motion  process.  We  begin  to  establish  this  claim  by  invoking 
Skorodhod's  theorem  (Skorodhod  (1965)  Section  6,  p.  9)  implying  in  this 
situation  that  (with  probability  one  in  some  probability  space)  we  may 
suppose  that  the  processes  Z^(w)  converge  uniformly  on  compact  w 

2 

intervals  to  a Brownian  motion  Z^w)  having  variance  parameter  <5  . We 
then  have 


lim  R.(s)  = lim{Z^(G^(s))  - Z^G^s))} 

X-*-°°  X-**> 

+ lim{Z  (G.(s))  - Z (s/k  LP  )} 
oo  x 00  eq 

X-x»  n 


+ lim{Z  (s/k  LP  ) - Z,(s/k  LP  )} 

°°  eC*  X e<5 


= 0, 


where  again  the  convergence  is  uniform  on  compact  s intervals  in  view  of 
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the  inequality  s/(tc  L^)  < Gx(s)  < s/(k  L^).  It  follows  that 

lim  Vx(s)  = -(1/K  Lgg)  ZJs/k  LPq) 

X-*°° 

= V (s) 

00 


where  V (s)  is  a Brownian  motion  process  having  variance  parameters 

oo 

Etvjl)2]  = S2/(k  LPq)3. 

2 

As  an  illustration,  we  evaluate  6 for  the  two  state  Markov  load 
process  L(t)  that  introduced  our  study  of  random  loads.  Recall  that  the 
possible  states  are  L + A,  and  denote  the  transition  function  for  L(t) 
by  Q in  the  form  Q+_(t)  = Pr{L(r)  = L - a|l(0)  = L + A), etc.  We  have 

Q ( a ) = Q++(a)  = i(l  + e"2Aa) 

and 

V(o)  = Q-+(o)  = l(1  ■ e_2XO)- 

Since  L(t)P  is  equally  likely  to  begin  in  either  of  its  states,  its  mean 

value  is  E[L(t)P]  = i(L+A)P  + ^(L-A)P  = LP  , and  for  a < t, 

i Z eq  ’ 

E[L(o)PL(t)P]  = (l/2K2)[a+a+Q++(T-o)  + a+a_  Q+_(t-a) 

+ a_a+  Q_+( r-o)  + a_a_  Q__(x-o)j 


t 
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which,  after  subtracting  the  product  of  means 
simplyfying  leads  to  the  covariance 


E[L(t )P  L(a)P]  = L2p  and 
eq 


E[F(L(ct))F(L(t))]  = k2  Cov[L (o)P,  L(t)P] 


b2  g-2X | t-o | 


where  3 = (a+-a  )/2  = C(L+A)P  - (L-A)P]/2. 

Then 

/J  /®  E[F(L(s))F(L(o))]do  ds  = g2  /J  /®  e'2A(s_0)  d0  ds 

= (B2/2X)  Jj  (l-e'2As)ds 
= (32/2X){T  + (2X)_1(e“2XT  - 1)} 

and 

62  = lim  | /T  fS  E[F(L(s))F(L(o))]do  ds = 32/X. 

That  is,  62  = X-1[(a+-a_ )/2]2  = X'^k  ((L+A)p  - (L-A)p)/2]. 

Returning  to  the  general  stationary  load  process,  let  us  consider 
large  bundles  (n-**>)  and  set  X = 0n  so  that  X becomes  infinite  with 
n.  We  have,  on  the  one  hand. 


/n  {S  -p}  = Z 
n n 


where  Z has  a limit  normal  distribution  with  mean  zero  and  variance  o ; 
n 


while  on  the  other  hand 


25 


V.  (s)  = /(0n){G  (s)  - s/(k  L^} ) 

A on  eq 


is  converging  weakly  to  a Brownian  motion  process  having  variance 


coefficient  52/(k  LPq)3.  Since  the  bundle  failure  time  Wn  is  given  by 


W = G„  (S  ),  we  write 
n 0n  n 


/n{Wn  - p/(k  Leq^ 


/n{G.  (S  ) - S /(<  LP  )}  + /n{S  -y}/(<  LP  ) 
0n  n n eq  n eq 


= U//9)  Vne(Sn)  ♦ Z„/<‘  L^) 


= <!//»)  Vne(„)  * ZJW  L^>  + tn 


where 


?„  ' ' V„6(u)}' 


Now  lim  S = ii,  while  V . is  converging  (uniformly  on  compact  intervals) 
n-M»  n n0 


to  the  Brownian  motion  V , and  Z is  converging  to  the  normally  distri- 


buted random  variable  Z.  Thus  £n  ■+■  0 in  probability  and  we  have  the 


weak  convergence 


^{W  - y/U  Lp  )}  =>  (1//T)  V (y)  + Z/(k  LP  ) . 
n eq  e4 


It  follows  that,  asymptotically  for  large  n,  the  standardized  failure 


time  /ii  {WR  - y/(<  LPq)}  is  normally  distributed  with  mean  zero  and 


variance 


( k L?J30  (<  L^q)2 


eq 
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We  see  that  random  loads  add  a term  to  the  variance  of  the  failure 
time  in  addition  to  decreasing  the  asymptotic  mean  failure  time  from 
y / ( k LP  ) to  y/(K  LPq). 

Let  us  suppose  that  nQ  is  the  number  of  fibers  needed  to  achieve  a 

prescribed  target  mean  cable  lifetime  under  a fixed  load  L,  and  that  n^ 

is  the  greater  number  of  fibers  needed  to  achieve  the  same  target  mean 

cable  life  under  a random  loading  L(i)  having  an  associated  L = 

{E[L(t  . As  shown  earlier,  we  have  the  equivalence  n.  /n„  = L /L 

eq 

so  that  the  ratio  L /L  provides  some  measure  of  the  increased  "cost" 

eq 

associated  with  random  loads.  This  ratio  is  sensitive  to  the  fiber  break- 
down rule  parameter  p and  the  marginal  distribution  of  the  load  process 

as  we  shall  see. 

2 2 

Let  A = E[(L(x)  - L)  ] be  the  variance  of  the  load  process.  When 
A is  very  small  relative  to  L and  P is  moderate,  so  that  pA 
remains  small,  we  have  the  Taylor  series  approximation 


L /L  = {E[(L(t)/L)P])1/P 
eq 

= {EC  ( 1 + )P])1/P 


^ r,  1 / . A. 2,1/P 

= {1  + ^p(p-l)(^)  } 


On  the  other  hand,  suppose  that  Lmax  is  the  essential  supremum  of  L(t) 

in  that  Pr{L(x)  > L - e}  > 0 for  any  e > 0 while  Pr{L(t)  <.  L } = 1. 

max  — max 

Then,  as  is  well  known  (see  e.g.  Taylor,  A.E.  (1958)  p.  91) 


lim  L = lim  {E[L( t)P]  }1/p  = L . 
p-x»  eq  p-*-°°  max 


1 


27 


That  is,  the  behavior  of  Lgq  for  large  values  of  p is  quite  distinct 

from  that  at  small  values.  To  obtain  a better  understanding  of  the  effect 

of  p on  the  ratio  L /L,  explicit  formulas  have  been  derived  in  several 

eq 

special  cases. 

(I)  Recall  that  when  L(t)  = L + i,  each  with  probability  1/2,  then 


L /L  = {i(l+A/L)p  + i(l-A/L)P}1/P. 
eq  2 2 


(II)  Suppose  that  L(x)  is  normally  distributed  with  mean  L and  variance 

2 

A . Then  for  p = 2m,  where  m = 1,2,...  we  have 


Leq/L 


■( 


e[  (i  + l(-!--)p]V/p 


£ Pf  (2Pk)E[(L~L~)2kj]1/P 

{ 1 * X **'>£*or  (S/I'>2k} 


2k  1 1/p 


(III)  Suppose  L(t)  is  uniformly  distributed  on  the  interval  [L-A,  L+A], 

2 2 

Then  E[L(t)]  = L and  Var[L(x)  = A = A /3,  while 


.1+P 


ErL(T)Pi  - _L  fL+A  xPdx  = (L+A)  - (L-A) 

ELL(x)  J - 2A  JL_A  x dx  2A(l+p) 


1+P 


Then 


WL  * ^2a(T+p7  C(1+a/l)1+P  - (1-A/L)1+P]}1/P 


2/3A(l  +P  ) 


[(1  + v^A/L)1+P  - (1  - /3A/L)1+P^  1/P. 
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L /L  as  a function  of  the  coefficient  of  variation  A/L  for  various 
eq 

values  of  p and  distributions  I,  II  and  III  is  depicted  in  Figure  3.  It 
is  readily  apparent  that  the  ratio  depends  strongly  on  the  form  of  the 
distribution  of  the  random  loads  and  cannot,  for  example,  be  adequately 
summarized  by  the  mean  and  variance  of  this  distribution  alone,  even  at 
modest  values  of  A/L.  This  dependence  on  the  load  distribution  becomes 
more'  acute  at  the  higher  values  of  p. 

It  is  clear  that  parallel  systems  whose  elements  exhibit  fatigue  of 
the  power  law  breakdown  type  should  be  designed  with  caution  when  randomly 
varying  loads  may  be  encountered. 


1 
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The  effect  on  cable  reliability  of  random  cyclic  loading  such  as  that 
generated  by  the  wave  induced  rocking  of  ocean  vessels  deploying  these  cables 
is  examined.  A simple  model  yielding  explicit  formulas  is  first  explored.  In 
this  model,  the  failure  time  of  a single  element  under  a constant  load  is 
assumed  to  be  exponentially  distributed,  and  the  random  loadings  are  a two 
state  stationary  Markov  process.  The  effect  of  load  on  failure  time  is  assumed 
to  follow  a power  law  breakdown  rule.  In  this  setting,  explicit  results 
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concerning  the  distribution  of  bundle  or  cable  failure  time,  and  especially 
the  mean  failure  time,  are  obtained.  Where  the  fluctuations  in  load  are 
frequent  relative  to  cable  life,  such  as  may  occur  in  long-lived  cables, 
it  is  shown  that  randomness  in  load  tends  to  decrease  mean  cable  life,  but 
it  is  suggested  that  the  reduction  in  mean  life  often  can  be  restored  by 
modestly  reducing  the  base  load  on  the  structure  or  by  modestly  increasing 
the  number  of  elements  in  the  cable . 


In  later  pages  this  simple  model  is  ascended  to  cover  a broader  range 
of  materials  and  random  loadings.  Asymptotic  distributions  and  mean  failure 
times  are  given  for  cable  elements  that  follow  a Weibull  distribution  of 
failure  time  under  constant  load,  and  loads  that  are  general  nonnegative 
stationary  processes  subject  only  to  some  mild  condition  of  asymptotic 
independence . When  the  power  law  breakdown  exponent  is  large , the  mean 
time  to  cable  failure  depends  heavily  on  the  exact  form  of  the  marginal 
probability  distribution  for  the  random  load  process  and  cannot  be  summarized 
by  the  first  two  moments  of  this  distribution  alone. 
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